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ABSTRACT 

We report the discovery of kinematic shock signatures associated with a locahzed radio jet interaction 
in the merging Seyfert galaxy NGC 5929. We explore the velocity-dependent ionization structure of 
the gas and find that low ionization gas at the interaction site is significantly more disturbed than high 
ionization gas, which we attribute to a local enhancement of shock ionization due to the influence of the 
jet. The characteristic width of the broad low-ionization emission is consistent with shock velocities 
predicted from the ionization conditions of the gas. We interpret the relative prominence of shocks 
to the high density of gas in nuclear environment of the galaxy and place some constraints of their 
importance as feedback mechanisms in Seyferts. 
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1. INTRODUCTION 

Most Seyfert AGN are associated with weak nuclear 
radio sources which show radio spectral indices and 
morphologies (if resolved) consistent with synchrotron- 
emitting jets and lobes (Ulvestad & Wilson 1984; Nagar 
et al. 1999), though their physical properties are poorly 
constrained. However, they are effective avenues of ki- 
netic and thermal feedback from the active nucleus and 
may play an important role in determining the evolution 
of the central spheroid. 

In the small fraction of Seyferts with kpc-scale ra- 
dio jets, several studies have uncovered clear signatures 
of interactions between the jet and surrounding gas, in 
the form of disturbed emission line profiles in the inner 
Narrow-Line Region (NLR) , as well as close associations 
between resolved jet structures and NLR gas (e.g. Whit- 
tle et al. 1988; Capetti et al. 1996; Falcke et al. 1998; 
Cooke et al. 2000; Cecil et al. 2002; Whittle & Wilson 
2004). Depending on the physical make-up of the jet, 
relativistic or ram pressure can drive fast shocks, com- 
pressing, sweeping up and altering the appearance of the 
NLR. Postshock gas, with a temperature of several 10^ 
K, is a source of ionizing photons which couple the shock 
properties to the ionization state of the surrounding gas 
(e.g. Dopita & Sutherland 1996). Shocks and winds can 
also change the distribution of ISM phases in the NLR 
by destroying and ablating clouds (Fragile et al. 2005). 
While the active nucleus is usually the dominant source 
of ionization even in strongly jetted Seyferts (e.g. Whittle 
et al. 2005), widespread shocks driven by a jet can alter 
the ionization of the NLR by affecting the properties of 
the ISM. 

Studies of shock structure and energetics (Dopita & 
Sutherland 1996; Allen et al. 2008) predict strong differ- 
ences between the emission line spectrum of dense post- 
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shock gas and gas that is ionized by, but not in direct dy- 
namical contact with, the shock (the precursor). There- 
fore, a clear signature of shock ionized gas is a difference 
between the line profiles of low and high ionization lines, 
which are preferentially produced by post-shock and pre- 
cursor gas, respectively (Whittle et al. 2005). 

In this work, we present an HST/STIS spectroscopic 
study of NGC 5929, a local Seyfert galaxy with a well- 
studied bi-polar radio jet. Previous ground-based spec- 
troscopic studies find evidence of a localized interaction 
between the jet and the near-nuclear emission line gas 
(Whittle et al. 1986; Wakamatsu & Nishida 1988; Taylor 
et al. 1989; Wilson & Keel 1989; Ferruit et al. 1997). The 
datasets and analysis methods used are briefly reviewed 
in §2. Direct shock features and a picture of the interac- 
tion are developed in §3 and §4. We discuss the role of 
shocks in Seyferts and AGN feedback in §5. NGC 5929 
has a systemic heliocentric velocity of cz = 2492 km 
s~^ based on the stellar aborption line redshift of Nelson 
& Whittle (1995), which corresponds to 161 pc arcsec"^ 
(Ho = 75 km s'^ Mpc'^). 

2. OBSERVATIONS AND REDUCTIONS 

2.1. STIS Dataset 

Our principal dataset is two medium-dispersion STIS 
G430M long-slit spectra covering the lines of the 
[O III]AA5007,4959 doublet and H/3 (HST Program: GO 
8253; PI : Whittle). The data was spectrophotometri- 
cally calibrated using the CALSTIS pipeline. For each 
spatial increment, the continuum was modeled using a 
lower order polynomial and subtracted. These spectra 
were combined to generate a continuum-free 2-D spec- 
trum from both STIS slits The equivalent width of H/3 
is always greater than 10 A in the NLR of NGC 5929, 
and ~ 140 A at the location of the interaction. Therefore, 
corrections to the emission line strengths and proflles due 
to Balmer absorption are negligible. 

2.2. HST and Radio Images 

Emission line maps of NGC 5929 in 
[OIIIJAA4959, 5007+H/3 were prepared from archival 
HST narrow/medium-band WF/PC-1 images (Program: 
GO 3724; PI: Wilson). Details of these maps can be 
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PC1 [0III]+H^ 
MERLIN 6 cm Contours 
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TABLE 1 

STIS OBSERVATIONS ^ 



Aperture 


Dataset 


Exp. (s) 


Nuc. Oflfset (") 


NGC 5929 A 


O5G403010 


1524 


0.198 


NGC 5929 B 


O5G403020 


600 


0.390 



Fig. 1. — Panels of images and spectra of the nuclear region of NGC 5929. The 2D spectra in Panels 2 and 3 correspond to STIS long-slit 
apertures A and B respectively, as indicated on the images plotted in Panels 1 and 4. The radio map (gaussian smoothed by 0''3 to bring 
out its structure) is plotted in Panel 1 in contours. 

surface brightness emission joins the hotspots to the ra- 
dio core along PA 61°, which is also the rough orientation 
of the elongated emission line distribution. 

The diverse imaging and spectroscopic datasets are 
showcased together in Figure 1. In Panel 1, the 
[O III]-|-H/3 emission line image is displayed in greyscale, 
with the STIS slit positions and radio map contours over- 
plotted. Three emission line regions, identified by Bower 
et al. (1994) are labelled A, B and C. The SW radio 
hotspot coincides with the edge of Region A, while the 
NW radio hotspot lies near Region B. From the broad- 
band image in Panel 4, the emission line features are re- 
vealed as parts of a possibly contiguous region, crossed by 
a dust lane that obscures the line and continuum emis- 
sion and modulates the appearance of the NLR. This 
dust structure is part of a filamentary network extend- 
ing to several kpc (Malkan et al. 1998). 

The [OIIIJA5007 line from the two medium-resolution 
STIS spectra are plotted in Panels 2 & 3 of Figure 1. 
The slits are oriented NE-SW from the bottom to the 
top of the panel, and a dashed horizontal line marks the 
reference position on the slit, i.e, the point along the slit 
closest to the nucleus of the galaxy. Solid vertical lines 
indicate the systemic velocity. 

The [O III] line from both STIS spectra shows clear 
velocity structure. A detailed treatment of the kinemat- 
ics of the NLR is not the aim of this Letter and will 
be addressed in a later paper. A brief description will 
suffice here. Slit A traverses the brighter emission line 
regions and therefore gives the best picture of the ionized 
gas kinematics. The velocity of the line peak increases 
along this slit from almost systemic at the reference po- 
sition to a maximum value of +185 km s~"'^at a nuclear 
radius of ~ 0'.'65 SW (104 pc). This broad gradient is 
mirrored, though less clearly, to the NE in both Slits A 
and B. The gradient between the nucleus and Regions B 
and C are similar in both slits, implying that both Re- 
gions are part of the same gaseous complex, bisected in 
projection by the dust lane. The FWHM of the [O Illj 
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found in Bower et al. (1994). In addition, a high S/N 
F606W WFPC2 archival image of the galaxy was used 
to trace the stellar and dust geometry in and around 
the NLR. A reduced and calibrated 5 GHz radio map 
of NGC 5929 was obtained from the MERLIN archive 
(http://www.merlin.ac.uk/archive/). This data was 
first presented in Su et al. (1996). 

The various HST datasets were registered to a common 
astrometric frame using image cross-correlation tech- 
niques. The dust lane in NGC 5929 crosses the NLR, 
effectively obscuring the true position of the nucleus. 
Lattanzi et al. (1997) use the astrometric information 
of stars in the WF/PC-1 images of the galaxy to com- 
pare the locations of the continuum peak and the unre- 
solved core in a MERLIN radio map. Assuming the radio 
core corresponds to the true nucleus, they find that the 
WF/PC-1 continuum peak is offset by O'.'l ± 0'.'05 NW 
from the radio core. We have included this small correc- 
tion in our final registered images (see Panel 4 in Figure 
1 for a sense of the amount of correction involved). 

3. THE NLR OF NGC 5929 
3.1. Descriptive Framework 

NGC 5929 is part of a major galaxy merger with 
NGC 5930 (Arp 90). Clear signs of tidal tails and 
disturbed gaseous kinematics are seen in ground-based 
two-dimensional spectra of the merger system (Lewis & 
Bowen 1993), and twisting, filamentary dust structure is 
visible in the nuclear region (Panel 4 of Fig. 1) 

The radio source has been imaged at high resolution 
with the VLA (Ulvestad & Wilson 1984; Wilson & Keel 
1989) and MERLIN (Su et al. 1996). It exhibits a 
triple structure, with two bright compact steep-spectrum 
hotspots on either side of a faint flat-spectrum core that 
is unresolved even in the highest resolution maps. Low 



line goes through a rapid transition from ~ 125 km s~ 
to greater than 200 km s~^ at a nuclear distance of 0'.'5 
in the brightest portion of Region A. Regions B and C, 
on the other hand, exhibit uniformly narrow line proflles. 

3.2. Ionization Conditions 
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Fig. 2. — The [OIII] image with MERLIN 6 cm contours over- 
layed (left) and the [OIII]/H/3 log ratio map (right). The location 
of the slit aperture is indicated on the image as a rectangle, corre- 
sponding to the spatial range of the ratio map. A grayscale lookup 
table is plotted in the bar above the map. A fiducial value of 0.43 is 
used to mask regions with no significant line emission, which helps 
to improve the visibility of the ratio map. 
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Fig. 3. — A comparison of the line profiles of the [O III] line 
(black) and line (red) in four parts of the NGC 5929 NLR. The 
lines have been scaled to a common peak value. The jet interaction, 
in the outer part of Region A, is associated with broad H/3, unlike 
any other part of the emission line region. 

In Figure 2, a map of the [OIII]A5007/H;9 Hne ra- 
tio from the Sht A spectrum is plotted, alongside the 
WF/PC-1 [O III] image, overlayed with the contours of 
the 5 GHz radio map and the boundaries of the STIS 
aperture. The ratio image was created by first rebinning 
the two-dimension spectrum of each line onto a common 
velocity range and then dividing them to generate a map 
of the [O III]/H/3 ratio as a function of velocity and slit 
position. This ratio of is quite sensitive to ionization 
state, yet insensitive to dust reddening. 

The [O III]A5007/H/3 ratio is in the range of 2 - 6, in- 
termediate between Seyferts and LINERs (Kewley et al. 
2006). Variations in ionization are seen across the NLR. 
The average excitation of Region B is towards the low end 
of the observed range, with average [O III] A5007/II/3 3, 
while the inner part of Region A is more highly ionized 
([O III] A5007/II/3 ^ 5). However, surrounding the po- 
sition of the bright knot in Region A, a sharp drop in 
the ratio is seen. Interestingly, this change in ionization 



is velocity dependent. In the high velocity wings of the 
lines, the [O HI] /H/3 drops to almost unity implying a 
very low ionization state, while the central core of the 
lines remain at a more modest ionization. This trend 
is different from most Seyferts where the high velocity 
gas tends to be of equal or higher ionization than low 
velocity material (e.g., Pelat et al. 1981). 

This difference between the average ionization of the 
NLR and Region A is brought out well by comparing the 
velocity profiles of [O III]A5007 and H/3 in different parts 
of the NLR, as is done in Figure 3. The width of the H/3 
line compared to [O III] is substantially higher in Region 
A, compared to both the inner part of the same Region 
and the profiles in Region B and C. 

What is the reason for this peculiar ionization be- 
haviour? We believe this results from strong shocks 
driven by the radio plasma into the ISM. The compressed 
gas in the post-shock cooling zone is expected to share 
the high velocities of the shock, yet is predicted to have 
[O HI] /H/3 in the range of 0.7 — 4.0, depending on the 
shock speed and local magnetic field (Dopita & Suther- 
land 1996). The enhancement of line emission from 
shocks and possibly a precursor region, over and above 
the emission from gas photoionized by the AGN, could 
then account for the appearance of the bright emission- 
line knot in Region A associated with the radio hotspot. 
A moderate contribution to the line emission from shocks 
is also consistent with the HST spectroscopic study of 
Ferruit et al. (1999). From the relative strengths of the 
[SII] AA4069, 4077 and AA6717, 6731 hues, which are pro- 
duced profusely in the post-shock cooling zone, they esti- 
mate high temperatures (greater than 20, 000 — 50, 000 K) 
and somewhat low electron densities (around 300 cm~^) 
for the [SII]-emitting gas. These values are quite reason- 
able for normal diffuse ISM that has been compressed 
by a shock front. Based on a simple double gaussian de- 
composition of the H/3 line integrated across the knot, 
we estimate an approximate linewidth of 420 km s~^ for 
the broad component, which includes the high velocity 
wings. This nicely matches the predicted shock velocities 
from Ferruit et al. (1999) which were based on compar- 
isons to shock ionization models of Dopita & Sutherland 
(1996). While it is difficult to directly associate shock 
speeds with integrated line kinematics, this broad con- 
sistency adds credence to a shock-driven origin for the 
emission of the knot in Region A. 

4. THE NATURE OF THE JET INTERACTION 

Our combination of high spatial resolution optical 
spectroscopy and imaging provide the best existing view 
of the compact jet-ISM interaction in NGC 5929. High 
velocity gas is only seen at the location of the south- 
western radio hotspot, while in other parts of the NLR, 
linewidths are narrow and a velocity gradient expected 
from virial motion is seen, consistent with ground-based 
slit and IFU spectroscopy (Keel 1985; Stoklasova et al. 
2009). 

We adopt a simple model for the jet interaction, fol- 
lowing that of Whittle et al. (1986). The radio jet plows 
through the ISM as it propagates outward, driving strong 
shocks at its head (the radio hotspot). As this shocked 
gas cools and gets denser, it becomes the [S II]-emitting 
gas described by Ferruit et al. (1999) and creates the 
broad H/3 emission. Allen et al. (2008) model the prop- 
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erties of post-shock g function of shock velocity 

Vsh, pre-shock density and magnetic field strength. Tak- 
ing a magnetized shock as fiducial, the [S II]-derived 
densities and Vsh ~ 400 km s~^(from the linewidth of 
the broad component) imply a pre-shock density of 
n ~ 30 cm~^ and a shock cooling length around tens 
of parsccs. The pre-shock density estimate has a large 
uncertainty and could be as low as 1 cm^^ if the shocks 
are weakly magentized, which would also decrease the 
cooling length significantly. However, since Seyfert jets 
are generally associated with mG level magnetic field, we 
adopt our estimates for magnetic shocks as most likely. 
The cooling length we derive is resolved in our HST im- 
ages and may explain the small offset of 0'.'15 between 
the location of the radio hotspot and the peak of the line 
emission in Region A. This interpretation is very tenta- 
tive, since patchy dust obscuration is the region may also 
cause such an offset. 

Su et al. (1996) estimate fiux densities, radio spec- 
tral indices and source sizes for the various compo- 
nents of the radio source, from which we calculate radio 
source cquipartition pressures, using the relations of Mi- 
ley (1980) and assuming that the radio emitting plasma 
has a filling factor of unity and an ion fraction a = 10, 
and that the radio spectrum extends from 0.01 GHz to 
100 GHz with a constant spectral index of -0.82. The 
cquipartition synchrotron pressure of the radio hotspot 
is 10^^ dyne cm~^, which may be compared to the ram 
pressure of the shock front, rUp n ~ 8 x 10~^ dyne 
cm~^. The two pressures are comparable, consistent with 
the view that the shocked gas surrounds and confines the 
radio hotspot. It is worthwhile to note that such pres- 
sures are considerably higher than those typically found 
in Seyfert radio sources, which may indicate that this is 
a relatively young jet ejection event. Eventually, inter- 
actions with the surrounding gas will confine the jet fiow 
into a static lobe, like those found in late-stage interac- 
tions (Capetti et al. 1996; Whittle & Wilson 2004). 

5. DISCUSSION; SHOCK SIGNATURES IN SEYFERTS 

Why do we see such obvious signatures of shocks in 
this object and not in others? This may be due to the 
weakness of the AGN: NGC 5929 has an absorption- 
corrected hard X-ray luminosity of 1.8 x lO''^ erg 
- low compared to average local Seyferts (Cardamone et 
al. 2007). Or perhaps the nuclear environment of the 
galaxy is dense and gas-rich due to its ongoing merger. 
This will lower the ionization parameter of nuclear radia- 



tion and produce a more compact emission line region. In 
both scenarios, the influence of the nucleus at larger radii 
will be relatively unimportant, making shock ionized line 
emission visible against the general background of cen- 
trally ionized gas. If this is indeed the case, it implies 
that radiative shock ionization from nuclear outflows is 
widespread in Seyfert NLRs, but is usually secondary 
to nuclear photoionization processes and only visible in 
low-luminosity Seyferts or those with dense nuclear en- 
vironments. A rich nuclear environment raises another 
possibility that the jet is impacting a dense molecular 
cloud, enhancing the shock luminosity. This may ex- 
plain why no clear shock signatures are seen around the 
NE radio hotspot, though the obscuration of the main 
dust lane prevents a direct view of this region. 

The luminosity of the broad H^ component is 8 x 10^^ 
erg s""'^. Following Dopita & Sutherland (1996), the H/3 
flux scales with the total radiative flux from the shock, 
with a weak dependence on Vgh, giving shock luminosities 
of 5 X 10^^ erg s~^ - comparable or slightly less than the 
total luminous output of the AGN (taking a X-ray bolo- 
metric correction of ^ 10). Using standard relationships 
from Osterbrock (1989), the H/3 luminosity can be used 
to derive the mass of ionized gas in the broad component: 
3 X 10'* Mq. If this mass of gas was accelerated to Vsh ■ it 
would have a total kinetic energy of 5 x 10^^ ergs. Taking 
the approximate acceleration timescale to be the crossing 
time of the region of size 0'.'3 (48 pc) at Vsh (around 10^ 
yr), a lower limit on the 'kinetic luminosity' of the jet is 
estimated to be 1.5 x 10''° erg s~^. This is few to several 
percent of the total AGN energy output. Given that jet 
outflows are a relatively common feature of Seyfert ac- 
tivity and couple strongly to the ISM through shocks, jet 
driven feedback can effectively carry as much as energy 
as the AGN photon luminosity to kpc scales, and transfer 
a tenth or less of this energy in the form of kinetic energy 
to the NLR. This can have important consequences for 
the suppression of bulge star-formation and the energy 
budget and dynamics of circum-nuclear gas. 
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